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ABSTRACT 

This  paper  presents  a  new  algorithm  for  both  broad¬ 
band  and  narrowband  adaptive  bcamforming  which  is 
robust  to  both  desired  signal  misalignment  and  array  imper¬ 
fections.  By  utilizing  available  a  priori  information,  the  new 
algorithm  minimizes  desired  signal  cancellation  by  introduc¬ 
ing  new  soft  constraints  which  are  incorporated  into  the 
LMS  adaptation  algorithm. 


1.  INTRODUCTION 

Adaptive  beamformers  can  be  used,  with  minimal 
a  priori  information,  to  greatly  improve  the  signal-to-noise 
ratio  of  a  signal  in  the  prescence  of  hostile  interference 
(jammers).  However,  when  the  actual  signal  angle  of  arrival 
does  not  coincide  with  the  estimated  angle  of  arrival,  or 
when  the  actual  geometry  of  the  array  differs  from  the 
geometry  assumed  by  the  array  designer,  then  the  perfor¬ 
mance  of  the  array  may  be  degraded  due  to  adaptation. 

This  phenomena  is  due  to  the  fact  that  if  there  is  misalign¬ 
ment  or  geometry  imperfections,  in  the  absence  of  a  pilot 
signal,  the  usual  adaptive  array  attempts  to  minimize  total 
output  power  regardless  of  the  type  of  signal  present.  While 
this  is  desirable  if  the  received  signal  is  mostly  jammer  and 
receiver  noise,  when  the  desired  signal  component  become.? 
large,  the  receiver  attempts  to  reject  it,  thus  significantly 
degrading  performance  [2,3]. 

The  most  common  way  of  avoiding  the  problem  of 
rejecting  the  signal  along  with  the  noise  is  to  impose  some 
sort  of  linear  constraints  on  the  ways  in  which  the  array  is 
allowed  to  adapt.  Many  methods  have  been  proposed  to 
implement  these  types  of  algorithms,  but  the  most  common 
are  the  Frost  algorithm  and  the  Generalized  Sidelobe  Can¬ 
celler  (GSC)  structure  [4, 5).  Both  of  these  types  of  arrays 
impose  constraints  on  the  array  patterns  which  they  may 
generate,  usually  either  fixed  direction  gain  constraints,  or 
constraints  on  the  derivative  of  gain  vs.  angle  at  a  fixed  look 
direction.  In  the  broadband  case,  these  constraints  may 
become  quite  complex  in  some  instances. 

The  array  processor  proposed  in  th'is  study  is  designed 
to  allow  the  cancellation  of  interferers  away  from  the 
demred  rignal,  while  being  very  robust  with  respect  to  both 
array  grometry  and  signal  direction.  This  is  achieved  by  use 
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of  the  available  a  priori  knowledge  of  the  signal  direction 
and  power  spectral  density  to  generate  soft  constraints  on 
the  processor  weights. 

2.  SYSTEM  MODEL 

For  the  remainder  of  the  paper,  the  structure  of  the 
adaptive  array  under  consideration  will  be  that  of  a  zeroth 
order  GSC.  This  structure  has  been  shown  to  converge  to 
the  same  steady  state  solution  as  the  Frost  constrained 
adaptation  algorithm  with  one  directional  constraint.  We 
will  assume  we  have  sources  emitting  only  plane  wave  sig¬ 
nals  (i.e.  far  field  approximation).  We  will  also  assume  one 
jammer  with  unknown  direction,  and  a  signal  whose  direc¬ 
tion  is  random,  but  known  a  priori  not  to  be  uniformly  dis¬ 
tributed  over  the  entire  (— 7r,?r)  angle.  In  addition,  we 
assume  that  we  have  some  knowledge  of  the  time  autocorre¬ 
lation  of  the  desired  signal.  It  will  later  be  shown  that  exact 
knowledge  of  these  quantities  is  not  required  for  good  per¬ 
formance.  We  will  also  assume  that  the  antenna  array  is 
linear,  with  uniform  X/2  spacing  between  antennas.  The 
gain  of  each  antenna  is  assumed  to  be  a  complex  Gaussian 
random  variable  with  inphase  and  quadrature  parts  i.i.d. 
with  a  distribution  N(l,<r,).  This  model  is  also  appropriate 
for  an  array  where  the  antenna  placements  are  random,  but 
close  to  linear. 

It  should  be  noted  that  the  assumptions  of  a  uniform 
linear  array  and  one  jammer  are  made  purely  for 
mathematical  expediency,  and  are  not  necessary  for  the 
results  of  this  paper  to  bold. 

The  sidelobe  cancellation  matrix  B  is  designed  to  be 
an  N  X  N-1  matrix  of  rank  N-1,  such  that 

BT-'O  (1) 

where is  the  N  dimentional  vector  of  delays  seen  for  a 
plane  wave  signal  coming  from  the  constndnt  direction. 
Normally,  the  array  is  assumed  to  be  pre— steered  so  that 
the  constraint  direction  is  broadside  to  the  array,  and  thus  : 

r-t  (2) 

We  must  define  a  new  received  vector,  taking  into  account 
the  delayed  versions  of  the  data.  With  this  in  mind,  we 
define  new  NK  dimentional  vectors  x»  and'sf  as: 
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Where  the  vectors  ‘5Z^i)  and  "s^i)  are  the  vectors  associated 
with  the  data  and  signal  delayed  by  i  sampling  periods, 
respectively.  We  must  now  defiM  a  new  beamforming  from 
the  matrix  B.  The  new  matrix,  B,  will  now  be  an  (N-I)K  x 
NK  block  diagonal  matrix,  of  the  form: 

B^diaglB,  ...,  B) 

,  The  output  of  the  GSC,  using  our  above  stated 

assumptions  on  the  gain  of  the  antenna  elements,  is  given  by 
the  expression: 

T‘  K/i  Gxt  -  w '  B  G7,  (3) 

where  \k/3  &  vector  with  N  ones  in  the  I</2  block  and 

zeros  everywhere  else,  and  ^  is  the  (N-1)I<.  dimensional 
vector  of  weights  used  in  weighting  the  outputs  of  the  N-1 
tapped  delay  lines  following  the  B  matrix. 

The  antenna  gains  will  be  incorporated  into  a  matrix 
G: 

G  is  the  (NK  x  NK)  block  diagonal  matrix  of  the  antenna 
guns,  and  can  be  written  as: 

G  *  diag(  G,  . . . ,  G ) 

and  A},-  is  a  N(0,o®|)  complex  random  variable.  It_should  be 
I  mentioned  that  a  given  realization  of  the  matrix  G  may  be 

I  dependent  on  the  directional  properties  of  the  data  vector 

Itf,  but  in  ail  cases  the  statistical  properties  of  the  matrix 
remain  identical. 

We  now  define  the  output  of  the  beamforming  matrix 
as  tf: 

P  Tf^BGl^i 

Tf  will  be  an  (N-1)K  dimentional  vector.  The  optimization 
problem  of  minimizing  the  total  output  power  of  the  GSC 
translates  to  the  problem: 

MIN;,  E[\Tk/2GX-^*^\^\  (4) 

^  The  solution  to  this  problem  is  the  well  known  Weiner 

Filter  solution.  The  optimum  weight  vector  ^  is  given  by 
the  equation: 

wf  (5) 

■  The  matrix  R,.  is  defined  as: 

And  the  vector is  of  the  form: 

“BGR„^G^  Ifc/j 

r  Where  R„  £  (1?(  Jtl  ^  ] 

And  ■j?!  ^  and  1?!  *  denote  complex  conjugate  transpose 
and  complex  copjugate  of  x^  respectively. 


3.  ROBUST  ADAPTATION  ALGORITHM 

If  the  signal  direction  vector Tis  not  exactly  equal  toT 
or  if  the  variance  of  the  antenna  gains  is  not  very  small,  the 
array  pattern  generated  by  (5)  will  attempt  to  null  out  a 
strong  desired  signal.  To  deal  with  the  problem  of  random 
antenna  gains,  Jablon  and  Widrow  [1]  have  proposed  using  a 
soft  constraint  type  of  optimization.  A  suitable  penalt)' 
function  is  chosen,  and  the  optimization  problem  is  then 
changed  to  minimizing  the  original  functional,  (  in  this  case 
^1  I  —  w'lT  P  ]  )  plus  some  penalty  constant 

times  the  penalty  function.  If  the  penalty  constant  is  chosen 
to  be  large  enough,  then  the  optimization  problem  will  have 
a  solution  very  close  to  that  of  a  hard  constraint  problem, 
where  the  constraint  is  defined  as  the  penalty  function  being 
forced  to  take  on  the  value  zero. 

We  wish  to  generate  a  penalty  function  which 
preserves  the  desired  signal  while  allowing  the  adaptive  pro¬ 
cessor  to  eliminate  interfering  sources.  In  a  previous  paper 
|6]  we  have  dealt  with  this  problem  in  the  narrowband 
situation.  We  wish  now  to  extend  this  technique  to  the  case 
of  broadband  arrays.  We  will  be  interested  specifically  in 
the  case  of  preserving  the  signal  in  the  presence  of  array 
imperfections  and  desired  signal  misalignment. 

We  will  define  our  penalty  function  p(^)  as  one  which 
is  linear  in  the  amount  of  power  that  is  transmitted  through 
the  cancellation  branch  of  the  GSC  due  to  the  signal  com¬ 
ponent.  The  leakage  of  the  signal  vector  into  the  cancella¬ 
tion  branch  is  given  by  the  expression: 

L(5^ )  "i^'BGs^  (6) 


We  will  lake  the  expected  value  of  the  power  of  the  leakage 
signal  as  our  penalty  function.  Thus,  p(^)  becomes: 

p(^)  =£[  w'B  (  GsTif^G^  )B^w  *  ]  (7) 


Bringing  the  ensemble  expectation  inside  and  utilizing 
the  assumptions  that  the  errors  in  the  antenna  gains  are 
independent  of  each  other,  and  of  the  possible  signal  direc¬ 
tions,  we  get  the  simplified  result: 

p(^)-^'B(o*,R'. -I-R^  )B’ir*  (8) 


Here  we  have  defined  two  new  matrices,  R'„  and  R„,. 
The  matrix  R'^  is  a  modified  time  autocorrelation  matrix 
of  the  signal.  It  is  a  block  NK  x  NK  matrix,  with  blocks  of 
size  N  X  N.  Each  block  is  an  identity  matrix  weighted  by  a 
time  autocorrelation  coenicient  of  the  ngnal.  The  ith 
matrices  off  the  diagonal  are  weighted  by  the  ith  time  auto¬ 
correlation  coefficient.  In  the  degenerate  case  of  a  nar¬ 
rowband  processor,!,  e.  no  tapped  delay  lines,  the  matrix 
R'„  becomes  simply  an  identity  matrix  multiplied  by  the 
desired  signal  power 

The  matrix  R^  can  be  thought  of  as  the  space-time 
autocorrelation  matrix  of  the  ^gnal  vector  S(.  The  matrix 
will  be  a  Hermitian  Toeplltz  block  Hermitian  Toeplitz 
matrix,  where  the  nth  off  diagonal  element  of  the  mth  off 
diagonal  block  is  given  by  the  expression: 


L'l  .-t  S;  •  . 


4.  RESULTS 


»/2 

f  R{fnT  +  nds\n{$))p(6)d0 

-»/2 

Where  R{ )  is  the  time  autocorrelation  function  of  the 
desired  signal,  p{ )  is  the  probability  distribution  of  the 
direction  of  arrival  of  the  signal,  T  is  the  unit  time  delay  in 
the  tapped  delay  lines,  and  d  is  the  time  for  a  wave  to  pro¬ 
pagate  between  adjacent  antennas. 

We  see  that  in  the  degenerate  case  of  no  tapped  delay 
lines,  the  integral  in  the  above  expression  becomes: 

»/2 

J  R(ndsin{0))p[0)d6 

-t/2 

And  in  the  narrowband  case,  the  autocorrelation  function 
/{( )  becomes  a  complex  exponential,  causing  the  integral  to 
take  a  form  similar  to  that  of  a  Fourier  transform. 

Utilizing  the  new  penalty  function,  the  new  soft  con¬ 
strained  optimization  problem  becomes: 

MIN;j  -2Re[^'BGR„GfT;^/2  1 

+  ^‘B  [R„ -t-c(  -f  R„,  )lB’w  *  (9) 

At  this  point,  we  note  that  the  term  is  quadratic  in  the 
vector  and  thus  any  minima  of  this  expression  will  be 
global.  The  matrices  R^,  and  R'^  are  assumed  to  be  known 
a  priori,  but  as  we  will  show  later,  any  reasonable  estimate 
will  suffice  for  the  proposed  algorithm. 

We  now  take  the  instantaneous  approximation  to  the 
gradient  vector  of  our  error  surface  using  the  usual  LMS 
methods  and  subtract  some  multiple  of  it  from  our  present 
solution  for  This  method  gi^  js  a  new  adaptive  algo¬ 
rithm  for  the  tap  weights; 

^*+1  -(I-2/i('if''u‘  +  -^B[o2,R'„+R«,  ]B’  ))wt 

+  2fStt  ‘Tx/2B'x,  '  (10) 

Or,  after  re-arranging  terms: 

-  (I-r?o*,BR'„B^  ->7BR„B’  [w* 

+  2#iy*'u'*  (11) 

Where  the  term  y*  is  the  output  of  the  array  at  the 
kth  iteration.  This  algorithm  can  be  shown  to  converge  for 
pomtive  values  of  ft  such  that; 

p,  < - — - 3—  (12) 

tr  I  R..  +  {  o*,R'«  -b  R«  )BM 

The  new  algorithm  requires  that  a  matrix  vector  multi¬ 
plication  take  place  during  each  tap  weight  vector  update. 
In  general,  this  would  be  an  order  ((N— 1)*K*)  operation. 
However,  since  the  matrix  in  the  multiplication  is  Hermitian 
ToepUtz  block  Hermitian  Toeplitz,  the  multiplication  can  be 
seen  as  a  convolution,  and  thus  it  can  be  computed  by  FFT 
in  order  ((N— l)Klog((N— 1)K)),  or  simplifying,  order 
((N— l)Klog(N— 1 )  -f  (N— 1  )K  log(  K )).  By  utilizing  this  fact, 
the  large  size  of  the  matrix  needed  for  the  robust  algorithm 
does  not  preclude  it’s  implementation. 


The  new  algorithm  requires  the  knowledge  of  the  signal 
space-time  autocorrelation  matrix  R^  and  the  desired  sig¬ 
nal  time  autocorrelation  function  R{)  a  priori.  In  practice, 
the  time  autocorrelation  function  may  be  well  known.  If  it  is 
not.  then  a  robust  estimate  of  the  autocorrelation  function 
may  be  used.  The  desired  signal  space-time  autocorrelation 
matrix  may  then  be  approximated  by  assuming  that  the  sig¬ 
nal  will  definitely  fall  within  the  main  beam  of  the  array  (or 
else  the  array  would  be  missteered),  and  then  taking  the 
maximum  variance  distribution  within  this  arc.  Using  this 
method,  the  spatial  distribution  of  the  desired  signal  is 
assumed  constant  within  some  arc  symmetric  about  broad¬ 
side,  and  zero  elsewhere. 

The  quiescent  pattern  (ignoring  misadjustment)  gen¬ 
erated  by  the  new  adaptation  algorithm  was  calculated 
analytically  for  a  broadband  array  with  11  antennas  and 
tapped  delay  lines  of  length  5  for  signals  having  a  PSD  given 
by  figure  2.  The  output  SINR  (signal  to  interference  plus 
noise  ratio)  was  calculated  as  a  function  of  jammer  angle  of 
arrival  for  a  fixed  signal  location,  and  again  as  a  function  of 
signal  angle  of  arrival  for  fixed  jammer  ■  location.  The 
results  show  that  even  with  inexact  knowledge  of  signal 
time  autocorrelation  and  direction  of  arrival,  the  new  robust 
algorithm  performs  quite  well,  almost  eliminating  the  effects 
of  signal  misalignment. 

We  then  simulated  both  the  robust  algorithm  and  a  Otli 
order  GSC  on  an  ll-by-5  processor  using  gaussian  data  hav¬ 
ing  identical  PSD  to  that  in  the  SINR  calculation.  The 
results  of  the  simulation  show  a  marked  improvement  of  the 
robust  algorithm  over  the  0th  order  GSC. 


5.  CONCLUSION 

We  have  developed  a  new  algorithm  for  robust  beam¬ 
forming  for  both  narrowband  and  broadband  arrays  in  the 
prescence  of  uncertain  signal  directions  and  array 
geometries.  The  new  algorithm  incorporates  available 
a  priori  information  into  a  penalty  function  which  imposes 
soR  constraints  on  the  adaptation  of  the  array.  The  new  soft 
constrained  algorithm  is  shown  to  perform  well  even  with 
reasonably  poor  information  about  the  derired  signal. 
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Fig.  1:  Broadband  GSC 


Fig.  2:  Power  spectral  density  of  simulated  signal.  Jammer. 
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Fig  3;  Output  SINK  vs.  Jammer  angle,  no  misalignment  Fig.  4:  Output  SINK  vs.  Jammer  angle.  0.02  rad.  mlaallgunent 
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Fig  5;  Output  SINK  vs.  Signal  angle.  Jammer  at  0.65  rad. 


Fig.  G:  Array  patterns  for  0th  order  CSC,  Robust  processor. 
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